Effect of coupling between junctions on superconducting current and charge 
correlations in intrinsic Josephson junctions 
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Charge formations on superconducting layers and creation of the longitudinal plasma wave in 
the stack of intrinsic Josephson junctions change crucially the superconducting current through 
the stack. Investigation of the correlations of superconducting currents in neighboring Josephson 
junctions and the charge correlations in neighboring superconducting layers allows us to predict 
the additional features in the current -voltage characteristics. The charge autocorrelation functions 
clearly demonstrate the difference between harmonic and chaotic behavior in the breakpoint region. 
Use of the correlation functions gives us a powerful method for the analysis of the current -voltage 
characteristics of coupled Josephson junctions. 



The intrinsic Josephson junctions (IJJ) have a wide 
interest today due to the observed powerful coherent ra- 
diation from the stack of IJJiA The radiation is related to 
the region in the current-voltage characteristics (CVC) 
closed to the breakpoint region (BPR)i^i^ii The resis- 
tively and capacitively shunted junction (RCSJ) model 
and its different modifications are well known to describe 
the properties of single Josephson junctions, giving a 
clear picture of the role of quasiparticle and supercon- 
ducting currents in the formation of CVCi^'^'i In the case 
of a stack of IJJ the situation is cardinally different. The 
system of the coupled Josephson junctions has a multi- 
ple branch structure and it has additional characteristics: 
the breakpoint current, the transition current to another 
branch and the BPR width. The breakpoint features 
were predicted theoretically and observed experimentally 
recently.^ We demonstrated that the CVC of the stack 
exhibits a fine structure in the BPR.i2. The breakpoint 
manifests itself in the numerical simulations of the other 
authors as wellfiS 

The breakpoint current characterizes the resonance 
point, at which the longitudinal plasma wave (LPW) is 
created in stacks, with a given number and distribution 
of the rotating and oscillating IJJ. These notions should 
be taken into account to have a correct interpretation of 
the experimental results. The investigation of the cou- 
pled system of Josephson junctions with a small value 
of the coupling parameter (as in the case of capacitive 
coupling), allowed us to understand in a significant way, 
the influence of the coupling between junctions on phys- 
ical properties of the system. The capacitive coupling 
is realized in nanojunctions if the length of the junction 
is comparable to, or smaller than the Josephson pene- 
tration depth at zero external magnetic field. Coupling 
between intrinsic Josephson junctions leads to the in- 
teresting features which are absent in single Josephson 
junction. Still, the superconducting current in the cou- 
pled system of Josephson junctions with LPW has not 
been investigated in detail; this includes the role of the 
correlations of the superconducting currents in different 
junctions and the charge correlations on superconducting 
layers. 



Here we study the phase dynamics of an IJJ stack in 
high- Tc superconductor. The CVC of IJJ are numeri- 
cally calculated in the framework of capacitively coupled 
Josephson junctions model with diffusion current i^^i^^ 
We find that the behavior of the superconducting current 
in the coupled system of Josephson junctions is essen- 
tially different from that of a single Josephson junction. 
It is demonstrated that superconducting current in the 
stack of IJJ reflects the main features of the breakpoint 
region; in particular, the fine structure in the CVC. We 
study various correlation functions in the characteristics, 
and observe that the correlations amongst the charge on 
different superconducting layers and correlations in the 
superconducting currents of different junctions lead to 
the detailed features of CVC in the BPR, and also pro- 
vides additional information about the phase dynamics 
in the IJJ. 
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FIG. 1: (Color online) The CVC and time averaged supercon- 
ducting current < sin > with no noise in current {In — 0) , 
and with noise amplitude I„ = 10~*. The inset shows an 
enlarged view of the BPR. 
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To find the CVC of the stack with N IJJ, we solve a 
system of N dynamical equations for the gauge-invariant 
phase differences (pi{t) = (pLi+i{t) = di+i{t) - ei{t) - 

JI~^^ dzAz{z,t) between superconducting layers (5*- 
layers), where 9i is the phase of the order parameter in 
the S-layer and Az is the vector potential in the barrier. 
In our simulations we use a dimensionless time r = tujp, 
where iOp is the plasma frequency LOp = \pleT^hC ^ Ic 
is the critical current, and C is the capacitance. The 
voltage is presented in units of Vq = ?kjp/(2e), and the 
current in units of I^- The system of equations has a 
form d'^LpijdT^ = Aui\l — smLpi — j3d(pi' /dr] with 
matrix A given in Ref. 4, for periodic and nonperi- 
odic boundary conditions (BC). Using Maxwell's equa- 
tion div{eE) = p/sq, where e and Eq are relative dielec- 
tric and electric constants, we express the charge den- 
sity Qi (we call it just charge) in the S-layer I by the 
voltages Vi and V;+i in the neighboring insulating layers 
Qi = QoctiVi+i - Vi), where Qo = eeoVb/r|), and ro is 
the Debye screening length. Solution of the system of 
dynamical equations for the gauge-invariant phase dif- 
ferences between S-layers gives us the voltages Vi in all 
junctions in the stack, and allows us to investigate the 
time dependence of the charge on each S-layer. The time 
dependence of the charge consists of time and bias cur- 
rent variations. We solve the system of dynamical equa- 
tions for phase differences at fixed value of bias current 
/ in some time r domain (0, T^) with a time step (5t; we 
then change the bias current by the current step 61, and 
repeat the same procedure for the current I + SI in the 
new time interval {Tm, '2-Tm)- The values of the phase and 
its time derivative at the end of the first time interval, 
are used as the initial conditions for the second time in- 
terval and so on. In our simulations we use = 25000, 
5t = 0.05, 51 = 10~^ and the dimensionless time tr is 
recorded as t^. = t -\- T^^Iq — I)/ 51, where /q is an initial 
value of the bias current. In this paper the CVC and the 
time dependence of charge oscillations in the S-layers are 
simulated at a = 1, /3 = 0.2, using periodic BC. The de- 
tails concerning the model and the numerical procedures 
have been presented beforej^i^ii^ 

Fig. [1] shows the CVC (green and blue curves) and 
the time averaged superconducting current < sin^j > 
(brown and red curves ) without random noise in cur- 
rent, and with noise amplitude 10~* for stack with nine 
IJJ. Our simulations of the CVC without noise gives us 
the value of the return current as Ir/Ic = 0.2517, which 
coinciding with the value obtained from the RCSJ model. 
In fact, in this model the relation between the return and 
critical currents for /3 << 1 has the form Ir/Ic = 4/3/7r, 
so that at /3 = 0.2 we get Ir/Ic ^ 0.2546. The super- 
conducting current without noise demonstrates the stan- 
dard increase before transition to the zero voltage s^^Xe^ 
The noise in current helps create the LPW in the stack 
and influences the superconducting current; the partic- 
ular value of noise is not very important. The creation 
of the LPW in the stack of junctions changes the CVC 
drastically, leading to the breakpoint and breakpoint re- 
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FIG. 2: (Color online) (a) The correlation of superconducting 
currents in the stack with nine IJJ: the color curves plot the 
correlation functions Cj^+i sin(/3-, (r) sin(/3-,+i(r) > as a 
function of bias current for all j, the black curve shows the 
corresponding CVC. Label {2,3}, for instance, indicates the 
charge or current correlations between layers 2 and 3. (b) 
The enlarged region in the current interval (0.561, 0.564). 
Notations for current and charge correlation functions are also 
given in the inset. 



gion in CVC. Inset shows the enlarged part of the Fig. [T] 
in the breakpoint region, where arrows indicate the coin- 
cidence of the main features of CVC and the supercon- 
ducting current. The point B shows the point where a 
charge appears on the S-layers; point Be is a breakpoint 
on the CVC and it reflects the breakdown in the sharp 
increase of the charge's value.— If we take a sum of all 
equations — (1 — aV^)(/ — /3(/3/ — sin(^/) for N IJJ 
and then find its average in time, we obtain the equation 
V ~ %{-^~ < sirup >). This equation clearly shows why 
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our results; namely, the V-curve and the < simp >-curve, 
demonstrate the same features. 

To investigate the origin of the CVC features in 



the BPR, we study the correlations, C: 



of 



perconducting currents in the neighboring junctions 
j and j + 1: ^^,,+1 sin93j(T) sin(y5j+i(T) >= 



lim 



(^j(r)sin (pj+i{T)dT, 



{T^-Ti] 

where the brackets <> mean averaging over time. The 
Cj as functions of bias current I / Ic are presented in 
Fig. [5^ for j — 1,...,9. All curves practically coincide 
in the interval starting from the breakpoint B till point 
Be- Then we observe regular, but different, behavior 
for different j till point Ci, and then again Cjj_^_i is 
nearly the same for all j in the chaotic region. As we 
can see in Fig. [2^, eight of these functions come in pairs 

(CI.S' ^2,3)' (^5,61 C'f_2); (C|_7,C| i), (C|_g,Cf_2); and 
one function, C| 4 stands by itself. Let us note that 
because we consider a periodic BC, the number assigned 
to a junction is only a label. The black curve shows the 
outermost branch of the CVC of the stack with nine IJJ. 
We can see that the features of the correlation functions 
coincide with the features of CVC, i.e. they manifest 
themselves in the CVC curve. In Fig. [2}d the part of 
Fig. [2^ is plotted on an expanded scale. Arrows show 
the points where the features of Cj coincide with 
those of CVC. In the inset to this figure we clarify the 
formation of the charge and superconducting current 
correlation functions by the diagram. Study of ^j^^ 
allows us to find new features in CVC which were not 
noticed in the previous studiesj^ Particularly, in Fig. [2] 
we indicate the points D2, -D3, D^, and where the 
curves of the correlation functions diverge, and points 
C/i, /72, [/a, t/4, where they converge. 
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FIG. 3: (Color onIine)ProfiIe of time dependence of the charge 
in the S-layers in the stack with 9 IJJ at a — 1, P = 0.2 and 
with periodic BC. 

To understand the origin of these features, we study 
the time dependence of charge in the superconducting 
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FIG. 4: (Color online) The correlation of charge in the stack 
with nine IJJ: the color curves plot the correlation functions 
C^ii+i =< <3;(''")Qi+i(''")) > for all I as a function of bias 
current; the black curve shows the corresponding CVC. The 
inset shows the enlarged region near the point Be- 



layers. In Fig. [3] we show the profiles of the charge os- 
cillations in all layers for the stack of nine IJJ. The odd 
number of junctions in the stack at periodic BC leads to 
the case where the charge dynamics in one layer, qualita- 
tively differs from the others; namely, the layer 4 demon- 
strates a specific time dependence of the charge oscilla- 
tions. We will refer to this layer as the "specific layer" 
or the sp-layer for short. The value of the charge on the 
sp-layer is smaller than on the other layers up to point 
Y, and practically staying at zero, in the interval from 
point T till point Y. Such difference leads to the maximal 



value of C: 



3,4- 



The correlation function C| 4 is different 



from the other correlation functions, because both phase 
differences in it ((^34 and (^45) are related to this unique 
layer. Any other correlation function has a partner which 
includes the same type of junction, for example (C| 5 and 
C| 3). Fig. [lb shows the enlarged region of Fig. in the 
current interval (0.561,0.564). It clearly demonstrates 
that the features of the correlation functions reflect the 
features of CVC . 

Let us investigate the charge correlations in the 
neighboring layers, using C^^i+i =< Qi{t)Qi+i{t)) >= 

lim(^T^^T,)^oo (rJ-Ti) It"' Qi{T)Qi+i{T)dT ■ Here we 
should emphasize that the index / is a layer number. 
Fig. |4] presents the dependence of Cf;^]^(///c) for 
different I for the stack with nine IJJ dX a = 1, (3 — 0.2, 
with periodic BC. The negative sign of the Cf;_|_]^ is 
due to the 7r-mode of charge oscillations, so that the 
product of positive and negative charges gives us a 
negative sign for this function. Again, the correlation 
functions Cfij^^ come in pairs: there are four pairs 
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FIG. 5: (Color online) (a) The CVC of the stack with nine 
IJJ; (b, c, d, e) The time dependence of the autocorrela- 
tion function Cf in the different parts of the BPR at / = 
0.57080, 0.56570, 0.56330, and 0.55800, respectively. 

(^4,5' ^3,4)^ (^5,6' (Cg 7, 2), (Cfs' C'g^i); here, 

the correlation function that stands by itself is C| g for 
layers 8 and 9, which are the farthest from the sp-layer 
(layer 4). The correlation functions (Cfg and CI4) 
between T and Y are close to zero, because the charge 
on the the sp-layer is close to zero in this interval. The 
inset shows the enlarged region around point Be, where 
the charge on the sp-layer demonstrates its specific 
dynamics (see Fig. [3]). The remarkable fact is that 
all pairs of current correlation functions C'jj+i and 
charge correlation functions Cfi_^^ form a loop reflecting 
this specific dynamics around this point! It means 
that the correlation functions reflect the correlations 
of phase dynamics amongst all layers and all junctions 
in the stack, even if the layers or junctions are far 
from each other. This is a demonstration of LPW in 
the system in other language. An interesting feature 
is observed in the chaotic region: At transition to the 



chaotic behavior (point Ci), the values of all correlation 
functions approach each other. The chaotic region will 
be discussed in detail elsewhere. 

To distinguish the harmonic oscillations in the BPR 
from the chaotic behavior, we study the autocorre- 
lation of the charge on the S-layers by the auto- 
correlation function =< Qi{t)Qi{t — ri) >= 

lim(T„_Ti)^oo (rJ-Ti) Qi('^)Qii'^ ~ Ti)dT. The au- 
tocorrelation function allows finding the repeating pat- 
terns (periodic signals) which have been buried un- 
der noise, or identifying the missing fundamental fre- 
quency in a signal, implied by its harmonic content. 
In Fig. [5^ we show the CVC highlighted at / = 
0.57080,0.56570,0.56330, and 0.55800, where the time 
dependence of the autocorrelation function Cf 
Q5{''')Q5{t — Ti) > is investigated. Results are presented 
in Figs. [5jb)- [5l^e). The different character of the auto- 
correlation function reflects the phase dynamics features 
in these parts of the BPR. At / = 0.5580, inside the 
chaotic domain, we can clearly distinguish the periodic 
motion from the chaotic one, in that C§ decays to zero 
with increasing time ti . 

In summary, the phase dynamics of intrinsic Josephson 
junctions in the high- Tc superconductors is theoretically 
studied. We establish a correspondence between the fea- 
tures of current- voltage characteristics, and the super- 
conducting current in the breakpoint region. We investi- 
gated the superconducting current in the coupled system 
of Josephson junctions with LPW and clarified the role 
of the superconducting current correlations in different 
junctions, in the formation of the total CVC. We demon- 
strated that the correlations of the superconducting cur- 
rents in neighboring junctions and the correlations of the 
charge on superconducting layers manifest themselves as 
the features on the CVC, as a consequence of the phase 
dynamics in the breakpoint region. We showed that the 
correlation analysis is a powerful tool for the investiga- 
tion of the CVC of the intrinsic Josephson junctions. 
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